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ABSTRACT 

As part of a JCMT Legacy Survey of star formation in the Gould Belt, we present early 
science results for Taurus. CO J=3-2 maps have been secured along the north-west 
ridge and bowl, collectively known as L 1495, along with deep 13 CO and C 18 J=3-2 
maps in two sub-regions. With these data we search for molecular outflows, and use 
the distribution of flows, HH objects and shocked H2 line emission features, together 
with the population of young stars, protostellar cores and starless condensations to 
map star formation across this extensive region. In total 21 outflows are identified. It is 
clear that the bowl is more evolved than the ridge, harbouring a greater population of 
T Tauri stars and a more diffuse, more turbulent ambient medium. By comparison, the 
ridge contains a much younger, less widely distributed population of protostars which, 
in turn, is associated with a greater number of molecular outflows. We estimate the 
ratio of the numbers of prestellar to protostellar cores in L 1495 to be ~1.3-2.3, and 
of gravitationally unbound starless cores to (gravitationally bound) prestellar cores to 
be ~1. If we take previous estimates of the protostellar lifetime of ^5xl0 5 yrs, this 
indicates a prestellar lifetime of 9(±3)xl0 5 yrs. From the number of outflows we also 
crudely estimate the star formation efficiency in L 1495, finding it to be compatible 
with a canonical value of 10-15%. We note that molecular outflow-driving sources have 
redder near-IR colours than their HH jet-driving counterparts. We also find that the 
smaller, denser cores are associated with the more massive outflows, as one might 
expect if mass build-up in the flow increases with the collapse and contraction of the 
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1 INTRODUCTION 



The Gould Belt Survey (GBS: IWard- Thompson et all 
l2007al ) is a large, legacy programme currently underway at 
the James Clerk Maxwell Telescope (JCMT). The goal is 
to gather both heterodyne line observations and submil- 
limetre (submm) continuum images of nearby star form- 
ing regions. The survey targets the Gould Belt, which is 
a ring of O-type stars and molecular clouds with a radius 
of ~350 pc. It is inclined at ~20° to the Galactic Plane 
and is centred on a p oint 200 pc from the Sun (ascendin g 
node, fi = 275 - 295; iTorra, Fernandez fc Figuerad 120001 ). 
The Gould Belt includes well known regions such as Orion, 
Taurus- Auriga-Perseus, Serpens, Lupus, Ophiuchus and the 
Pipe nebula. Taurus is one of four regions initially targeted 
in 12 CO J=3-2 (hereafter CO 3 -2), 13 CO 3-2 and C ls O 3-2. 
These four regions - Orion A l|Buckle et alj|20ich . Serpens 
(|Graves et al. I20T0I ). Ophiuchus (White et al., in prep.) and 
Taurus (this paper) - were selected as being complementary, 
active, nearby, relatively well known, and well studied. Ex- 
tensive submm continuum mapping, at 450 /xm and 850 /xm, 
is also planned as part of the survey (|Ward-Thompson et al.l 
l2007al ). 

In studies of nearby, low-mass star forming regions CO 
3-2 emission samples the denser, warmer gas which often 
envelopes cores where new stars are forming (temperatures 
~ 10-50 K; gas densities of the order of 10 4 -10 5 cm -3 ). CO 
3-2 is also a proven tracer of outflow activity (e.g. Hatchell, 
Fuller & Ladd 1999; Knee & Sandell 2000; Davis et al. 2000a; 
Hatchell, Fuller & Richer 2007a; Lee et al. 2007; Bussmann 
et al. 2007; Yeh et al. 2008). With the JCMT, moderate spa- 
tial resolution is attainable at 345 GHz (~14"), which often 
allows one to disentangle multiple outflows in clustered re- 
gions. The higher critical density for excitation of the 3-2 
line (over, say, the 1-0 emission line) can also lead to better 
maps of the dense, more collimated flow components. Molec- 
ular outflows are known to emanate from the very youngest 
sources, objects that are still accreting most of their final 
mass. New, large-scale surveys in CO can therefore be used 
to identify the locations of these embedded protostars, lead- 
ing to a more complete census of the youngest stars, and 
an indication of the overall youth and star formation activ- 
ity in a region (e.g. . Hatchell ct al. 2007a; D avis et al.ll2008l 
lHatchell fe Dunham|[2009l ; iBuckle et al]|2010l ) . 

Taurus is a much-studied region of low mass star for- 
mation (see Kenyon, Gomez & Whitney 2009, for a detailed 
review). Th e Taurus molecular cloud covers an area in excess 
of 100 deg 2 ijUn gcrcchts & Thaddcus 1987; Narayanan et al.l 
|2008| ; IGoldsmith et al.ll2008l ). LDN 1495 (L 1495 hereafter) 
is one of three main areas associated with a high density of 
young stars. L 1495 lies in the north-west corner of Taurus 
and coincides with Barnard dark nebulae B 7, 10, 209, 211, 
213 and 216. In molecular line and extinction maps L 1495 
comprises an extended, knotty filament - here referred to 
as the "south-east ridge" - and, at its northern ext r emity , 
a more diffuse cloud - the "bowl" . IGoldsmith et all (|2008l ') 
refer to the ridge as B 231 and label only the bowl L 1495. 
From their extensive CO observations they estimate masses 
of 1095 Mq and 2626 Mq for the ridge and bowl, respec- 
tively, and assign areas of 13.7 pc 2 (2.3 sqr degrees) and 
31.7 pc 2 (5.3 sqr. degrees) for each region (see their Table 
4). They also note that, quite remarkably, the SE ridge is 75' 



(3 pc) long yet only 4.5' (0.2 pc) wide. A similar collimated 
feature in the N GC 6334 star forming region has recently 
been analysed by [Matthews et all (|2008l ). 

In this paper we split the bowl into four areas, 
L 1495 N, S, E and W, and label the ridge L 1495 SE. 
L 1495 E harbours t he classical T Taur i star (TTS) CoKu 
Tau-1 and HH 156 (|Eisloffel fc Mundtl Il998l ); L 1495 W 
(B 209) is associated with the TTS CW Tau (Elias 1, 
Hubble 4) w hich drives HH 220 and possibly also 
HH 826-828 (|G6mez de Castrd 1 19931 ; iMcGroartv fc Ravi 
|2004| ). HH 390/391/392 and HH 157 are located along the 
L 1495 SE ridge (Gomez, Whitney & Kenyon 1997; Eisloffel 
& Mundt 1998); HH 157 is driven by the TTS Haro 6-5B 
(FS Tau B) and comprises a spectacular HH jet and bow 
shock that extend over 30"-40". 

L 1495 W has been mapped in C 17 1-0, C ls O 1- 
0, CS 2-1, N 2 H+ 1-0 (at FCRAO) and in NH 3 (1,1) 
and (2,2) emission (at Effelsberg) by iTafalla et all (|2002l ) 
- the region they label "L 1495" roughly coincides with 
L 1495 W. N 2 H+ 1-0 map s of th e south-east ridge are pre- 
sented bv lTatematsu et all (|2004l ) . More exten sive 13 CO and 
C 18 1-0 maps of Taurus are discussed by iMizuno et al.l 
1 19951 ). and lOnishi et all (|l996l ); CO and 13 CO 1-0 maps, 
covering 98 square degrees with 45" resolut io n, ha ve re- 
cently been pu b lished by iNaravanan et all (|2008l ) and 
IGoldsmith et al] l|2008l ). Finally, extensive yet high spatial 
resolution H 13 CO + 1-0 observations of the molecular cores 
in Taurus, includ i ng the L 1495 ridge and bowl, are discussed 
by lOnishi et all |2002i . hereafter referred to as OMK02). 
These H 13 CO + observations, to which we refer throughout 
the paper, were obtained at the 45 m Nobeyama telescope 
wit h a beam width of 20" . 

iKenvon &: Hartmannl (| 19951) present a list of ~300 
Young Stellar Objects (YSOs) in Taurus, derived from 
multi-wavelength observations and compl emented recently 
with data from the Spitzer Space Telescope l|Hartmann et all 
2005; lLuhman" ct al. 2006, though note that the Spitzer ob- 
servations do not cover L 1495 W or the bottom of the 
south-east ridge). Many of the youngest sources have been 
observed in the 1.3 mm continuum survey obtained by 
iMotte fc Andrei (|200ll ). In L 1495 the young stars are clus- 
tered toward the eastern and western sides of the L 1495 
bowl, and are tightly bound along the narrow, south-east 
ridge. The pr otostellar population as a whole is discussed in 
the review of IKenvon et all {2009), who also present an up- 
dated list of YSOs in Taurus. Kenyon et al. note in their ar- 
ticle that classical TTSs typically have near-infrared (near- 
IR) colours H — K s < 1 — 1.5 and mid-infrared (mid-IR) 
colours [3.6]-[4.5]~0. 5-1.0 and [5.8]-[8.0]~0.25-0.75; embed- 
ded protostars (Class 0/1 sources) are expected to be redder. 

In this paper we use an up-to-date list of "Taurus Young 
Stars", kindly provided by K. Luhmann, to search for can- 
didate outflow sources in L 1495. This list i s esse ntially the 
same as that presented by IKenvon et all {2009). We also 
assume that "starless" cores are low mass (M < 10 Mq) 
dense cores without comp act luminous sources of any mass 
(|Di Francesco et alj |2007). "Prestellar" cores are a subset of 
starless cores, since they must also be gravitationally bound 
I Ward- Thompson et aill2007bl ). Establishing whether a core 
is bound or not is observationally very difficult (although the 
virial theorem gives some indication). We therefore refer to 
cores that are not obviously associated with a young star as 
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Figure 1. Left - an overview of L 1495 in the north-west corner of Taurus. The image sh ows integrated C O 1-0 emission from the 
FCRAO survey of Goldsmith et al . (2008). The positions of the H 13 CO+ cores identified bv lOnishi etaD l|2002l l are marked with circles 
and labelled using their core numbering scheme. Right - a complete map of all regions observed with HARP in CO 3-2 (intensity 
integrated between 2.0 and 11.0 km s — 1 ), with the positions of the H 13 CO+ cores again marked. The areas mapped in 13 CO and C ls O 
3-2 are also indicated; maps of 13 CO emission, integrated between 4.9 and 9.1 km s _1 (L 1495 E) and 4.5 and 8.2 km s~ 1 (L 1495 SE), 
are shown inset as grey-scale images. In the main text we refer to L 1495 N, E, S and W as the bowl, and L 1495 SE as the south-east 
ridge. 



Table 1. Log of CO 3-2 observations in L 1495. 



Region" Area" Cube" Date" RA 6 Dec 6 PA C Map size c v d 

(J2000.0) (J2000.0) (kms" 1 ) 



Bowl 


L 1495 E 


L 


1495 


2007-11-24 


4 


18:04.8 


28:21:45 


90° 


20'xlO' 


6.6-7.1 


n 


// 


L 


1495-S 


2007-12-22 


4 


17:50.0 


28:09:00 


0° 


10'xl5' 


5.9-6.7 


ii 


// 


L 


1495-S2 


2008-09-17 


4 


18:48.0 


28:14:30 


0° 


10'xlO' 


~5.9-6.7 


ii 


L 1495 N 


L 


1495-NE1 


2008-11-08 


4 


16:59.0 


28:30:13 


0° 


10'xl4' 


~6.7 


ii 


// 


L 


1495-NE2 


2009-01-30 


4 


16:28.1 


28:45:33 


0° 


10'xl4' 


~6.7 


ii 


L 1495 S 


L 


1495-S3 


2008-11-08 


4 


17:34.0 


27:49:20 


0° 


10'x22' 


6.6-6.8 


ii 


// 


L 


1495-SE1 


2008-03-29 


4 


18:09.1 


27:35:17 


0° 


5'x5' 


5.6 


ii 


L 1495 W 


L 


1495-W 


2008-08-21 


4 


14:10.0 


28:12:00 


0° 


10'xlO' 


6.7 


ii 


// 


L 


1495-W1 


2009-01-31 


4 


15:54.7 


28:16:33 


0° 


22'xlO' 


~6.7 


ii 


// 


L 


1495-W2 


2008-10-07 


4 


14:34.0 


28:04:00 


0° 


13' x 13' 


6.7 


ii 


// 


L 


1495-W3 


2009-01-30 


4 


13:35.0 


28:22:13 


0° 


13' x 13' 


6.7 


Ridge 


L 1495 SE 


L 


1495-SE2 


2008-03-15 


4 


18:41.3 


27:22:00 


147° 


8'xl5' 


5.3-5.7 


it 


// 


L 


1495-SE3 


2007-11-25 


4 


19:42.4 


27:13:24 


147° 


8'xl5' 


6.5 


ii 


// 


L 


1495-SE3b 


2008-01-12 


4 


20:10.0 


27:05:00 


147° 


10'xlO' 


6.4 


ii 


// 


L 


1495-SE3c 


2008-09-17 


4 


19:22.0 


27:04:06 


147° 


9'x23' 


~6.4 


ii 


// 


L 


1495-SE4 


2008-08-21 


4 


21:14.1 


27:01:11 


113° 


10'x20' 


6.6-6.7 


ii 


// 


L 


1495-SE4b 


2008-10-07 


4 


22:22.0 


26:55:00 


113° 


10'xlO' 


~6.6 


ii 


// 


L 


1495-SE5 


2008-11-08 


4 


23:44.5 


26:39:00 


113° 


9'x 25' 


6.6-6.7 



"The areas in the bowl and ridge regions of L 1495 are made up of multiple data cubes, often observed on different nights. 
''Coordinates of the cube centres. 

c Cube map size and position angle of the map long axis (measured E of N); with "basket- weave" mapping the scan direction is along 
and orthogonal to this axis. 

d The LSR velocities of the H 13 CO + cores found in or near (if an approximate value is given) each region by OMK02, which we adopt as the local 
ambient gas velocity. 
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starless; cores that do appear to contain a young star are 
referred to as "protostellar" . 

Traditi onally, a distance of 140 pc has been used 
for Taurus (|Ungerechts fc Thaddeusill987l; lElias! Il978l). o n 
the basis of estimates from star counts 1 McCuskevI 19391 ) . 
the reddening of field stars as a function of distance 
l|Gottlieb fe Upsonl [l969). and photometric distanc es mea- 
sured to the exciting stars of reflection nebulae (|Racinel 
1968). Analysis of the velocity field in wide- field molecu- 
lar line studies suggests that ambient gas velocities change 
by only a few km s" 1 ac ross Taurus, and are particu- 
larly limited acro ss L 1495 dlJngerechts fe Thaddeuj|l987l ; 
iNaravanan et ah l l2008l ; lGoldsmith et al.ll2008l . see also Sect. 
3.1 below). This supports a more-or-less common distance 
to the molecular clouds and star forming regions in Taurus. 
Recent trigonometric parallax measurements suggest that 
the eastern portion of Taurus may be furthest from us, at 
a distance of ~ 161 pc, with the south (around T Tau) at 
an intermediate distance of 147 pc and the west, L 1495, 
marking the near s ide of the cloud at a distance of 130 pc 
jTorres et al.ll2009l ). However, these distances are based on 
parallax measurements for only a handful of stars. We there- 
fore adopt the canonical distance of 140 p c for L 1495, con- 
sisten t with previous studies, most notably iGoldsmith et al.l 
l|2008h and OMK02. 

Overall, our goal with this paper is to search for outflows 
in L 1495, identify their driving sources, and compare the 
properties of these outflows with those of the driving sources 
and associated molecular cores. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Data Acquisition 

As part of the JCMT Gould Belt Legacy Survey, the 325- 
375 GHz heterodyne receiver array HARP was used to map 
18 near-contiguous regions across the L 1495 bowl and along 
the south-east ridge in CO 3-2 emission (Fig.[T]and Table[T]), 
as well as two smaller sub-regions in 13 CO 3-2 and C 18 3-2 
emission. HARP is a 16 receptor (4x4 array) single side- 
band SIS receiver system acting as a fro nt-end to the AC- 
SIS digital auto-correlation spectrome ter (|Hovev et al . 2000; 
ISmith et al1l2003l ; iBuckle et al.ll2009l ). 

CO 3-2 data at a rest frequency of 345.79599 GHz were 
acquired in dual sub-band mode. Maps with modest and 
high spectral resolution were acquired simultaneously, with 
1024 ACSIS channels covering a bandwidth of 1.0 GHz at 
a spectral resolution of 977 kHz (0.85 kms -1 ), and with 
4096 channels covering a bandwidth of 250 MHz at a spec- 
tral resolution of 61 kHz (0.05 kms" 1 ). Only the latter are 
presented in this paper. In separate observations, 13 CO 3- 
2 (330.58796 GHz) and C 18 3-2 (329.33055 GHz) data 
were acquired simultaneously in dual sub-band mode using 
only the high spectral resolution set-up. At these frequen- 
cies the telescope Half Power Beam Width measures 14", 
corresponding to 0.0095 pc at a distance of 140 pc. 

Maps were obtained by scanning the telescope along 
rows parallel to the sides of the pre-defined mapping areas, 
each row separated by half or a quarter of an array to im- 
prove the spatial sampling. This map is then repeated but 
by scanning in a perpendicular direction. This strategy of 



"basket weaving" helps to generate natter maps and evens 
out the noise in the final data cube. For most of the data pre- 
sented here, 2-4 receptors were inoperable; for those fields 
where 4 receptors were unavailable (regions S2, W2, SE3c 
and SE4b) the array was stepped by a quarter of the array 
to prevent empty spaces or poorly sampled regions in the 
final maps. 

Depending on map size, the 12 CO data cubes were ob- 
tained by integrating for 0.6 sec or 1.2 sec (the "sample 
time") for each 7.3" pixel on the sky. A reference position 
(usually observed at the end of each row in a scan), located 
at 4 h 12 m 24.9 3 (J2000.0) 26° 50' 39", was used throughout. A 
sample time of 1.0 sec was used with the smaller 13 CO/C ls O 
maps. 

The CO 3-2 data were observed in grade 3 weather 
(t(225 GHz) ~ 0.08-0.12); better weather (r(225 GHz) ~ 
0.05-0.065) was used for the 13 CO 3-2 and C 18 3-2 data. 
System temperatures for the 12 CO data ranged from 340 
to 580 K (<450 K in most regions); for the isotopologues 
temperatures were typically 300-500 K. 

The regions mapped in 12 CO are listed in Table Q] 
these were defined based on large-scale extinction maps 
|Dobashi et al.ll2005l ). the surveys of dense cores conducted 
by OMK02, and the locations of Herbig-Haro (HH) ob- 
jects. Only regions with bright 12 CO emission were mapped 
in 13 CO and C ls O: two overlapping maps centred near 
H 13 CO+ cores 6, 7 and 11 were repeated a number of times 
to reduce the noise in the final cube; the full map covers an 
area of ~22'xl2', although note that the noise is roughly a 
factor of two higher in the north and west portion (~60%) 
of the map. Similarly, two overlapping maps covering cores 
13a, 13b and 14 were repeated in the L 1495 SE ridge; in this 
case, the map covers a rectangular area of approximately 
20' x 6', though again the north-west quarter of the map is 
about 3-times noisier than the rest of the data cube. In all, 
eleven overlapping maps were co- added in L 1495 E, and 
eleven in L 1495 SE. The regions covered in 13 CO and C 18 
are marked in Fig. [T] 



2.2 Data processing 

The individual CO 3-2 data cubes listed in Table[T](third col- 
umn) were reduced usin g the automated ORAC-DR pipeline 
l|Cavanagh et alj|2008l ). Bad or overly noisy spectra were 
identified by the pipeline, based on quality assurance crite- 
ricQ stipulated for the Gould Belt Survey legacy data, and 
these were first removed and noisy regions at the band edges 
trimmed. The resulting time-series data were re-gridded to 
a 3D data cube and regions of line emission identified and 
masked out. A fifth-order polynomial was fit to the result- 
ing emission-free regions and subtracted from each reduced 
cube. Cubes observed at the same spatial and spectral po- 
sition were combined and a new baseline region mask was 
determined from this higher signal-to-noise cube. This was 
applied to the individual cubes, which were then mosaicked 
and co-added together using a light Gaussian smoothing 
kernel with a Full Width Half Maximum (FWHM) ~7.3", 



1 Receiver system temperature >600 K, noise across the spec- 
trum exceeds 50% of the mean noise level in all spectra, or spec- 
trum baseline deviates by more than 10% from the mean. 
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40 20 4:20:00 40 20 19:00 



L 1495 E (12C0 3-2) 

Figure 2. 13 CO 3-2 channel maps across L 1495 E, plotted be- 
tween 4.8 and 8.0 km s _1 , in 0.4 km s — 1 intervals. H 13 CO+ cores 
are marked and labelled bottom-left; the small cross in the centre 
of each map is drawn as a reference marker. The velocity of the 
channel is marked top-right in each panel. 



to produce the final cubes and images. The inverse of this 
baseline mask identifies the line emission regions and this 
was used to determine the integrated intensity and veloc- 
ity images presented in this paper. Note that the Gaussian 
smoothing results in a spatial resolution to point sources of 
about 16". 

Special care was taken to remove st riping from the data 
cubes when present l|Curtis et al.| [2010). This was done by 
determining the receptor-to-receptor response for the whole 
time-series, and dividing that response out. This technique 
does assume a priori that each receptor is exposed to similar 
emission during the scan. Given the diffuse and extended na- 
ture of the CO emission from Taurus, this technique worked 
well in eliminating or, at the very least, in minimising the 
striping. 

fn CO 3-2 the noise is relatively uniform across the 
extensive region mapped. In the un-binned data (at the 
full 0.05km s _1 resolution), at velocities within a few km 
s" 1 of the ambient, the root-mean-square (RMS) noise level 
{T\(rms)) measures 0.10-0.20 K in L 1495-E, N and S (be- 
ing highest in the south), 0.08-0.15 K in L 1495-W (highest 
in the north-west), and 0.08-0.16 K along the south-east 
ridge (it is a little higher, ~0.22 K, in the data around cores 
10, 10b and 12a). 

The 13 CO and C ls O data were reduced in the same 
way. As noted above, the maps in 13 CO and C ls O were re- 
peated a number of times. We found that in order to produce 
the optimal coadded map, the worst quality data had to be 
discarded before being included in the map reconstruction. 
This was done by measuring the RMS of each of the spectra 
in the input data and masking those with the highest RMS 
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L 1495 SE (1200 3-2) 

Figure 3. 13 CO 3-2 channel maps across L 1495 SE, plotted 
between 5.15 and 8.15 km s — 1 , in 0.5 km s _1 intervals. H 13 CO + 
cores are marked and labelled top-right; the small cross in the 
centre of each map is drawn as a reference marker. The velocity 
of the channel is marked top-left in each panel. 



before generating the map. An iterative method was used 
to determine the optimal RMS value to apply as a mask. 
In each step of the iteration, the data were masked using 
a different limiting RMS before being recons tructed into a 
map using the Starlink MAKECUBE routine |jenness et al.l 
2008). The RMS in the reconstructed map was then mea- 
sured. The limiting RMS was varied to produce a recon- 
structed map with the lowest RMS. This map was taken to 
be the optimal map. In the final, coadded maps the RMS 
noise level measures Tj(rms) ~0. 13-0. 20 K in 13 CO and 
~0. 16-0.25 K in C ls O at the unbinned spectral resolution of 
0.05 km s _1 . As mentioned earlier, the noise level is some- 
what higher in the north-west regions of each map. 

In this article all images are presented in units of an- 
tenna temperature (TjQ. When calculating outflow parame- 
ters, integrated antenna temperature has been converted to 
main beam brightness temperatu re, Tj, = T* A / '(r)mb) , using 
an aperture efficiency r\ m b = 0.63 l|Buckle et al.ll2009T ) . 



3 RESULTS 

3.1 Overall distribution of CO 3-2 emission and 
large-scale velocity structure 

The entire region mapped in CO 3-2 with HARP is shown 
in Fig.[T] where we also show integrated 13 CO 3-2 images of 
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Figure 4. A plot of intensity-weighted radial velocity (first mo- 
ment) showing the distribution of blue and red-shifted CO 3-2 
emission around L 1495 E and L 1495 S. Black/blue through 
to red/white display velocities in the range 6.0 to 8.0 km s — 1 
(green/yellow represent intermediate velocities). Circles, crosses 
and triangles mark the positions of H 13 CO + cores (labelled), 
YSOs and HH objects, respectively. Integrated CO 3-2 emission 
(zeroth moment, from Fig. [TJ is over-plotted as contours; levels 
measure 4,6,8,10,12 K km s -1 . Note that there is CO 3-2 line 
emission throughout this entire map. 



the L 1495-E and L 1495 SE regions observed in this isotopo- 
logue (and in C ls O - data not shown). For comparison pur- 
poses, we also present the 13 CO 1-0 map of lGoldsmith et aU 
(I200SI) . 

When mapping L 1495 in CO 3-2 we followed the high- 
extinction ridge of dense cores obse rved in C ls O 1-0 and 
H 13 CO + 1-0 by lOnishi et all l|l99rj) and OMK02, respec- 
tively; the ridge is evident in the 13 CO 1-0 integrated in- 
tensity map reproduced in Fig. [T] In L 1495 the CO 3-2 
emission is generally faint, diffuse and/or optically thick: 
the integrated antenna temperature is < 10.0 K over 99% of 
the mapped region, while we estimate the 12 CO opacity, T12, 
to be in the range 3-38 in regions where 13 CO has been ob- 
served, (although values will be overestimated where profiles 
are self-absorbed - see Section 3.2 for details). There are a 
number of compact knots superimposed on this diffuse emis- 



sion, most of which are associated with outflows rather than 
cores. Indeed, the CO 3-2 emission features are by-and-large 
unrelated to the compact cores identified in H 13 CO+. Unfor- 
tunately there is very little submm continuum data in L 1495 
for comparison purposes (in their pointed 1.3 mm survey of 
young stars, Motte & Andre [2001] observe only nine sources 
in L 1495); only observations of cores 13a, 13b and 14 were 
found in the SCUBA (Submill imetre Common User Bol ome- 
ter Array) Legacy Catalogue (|Di Francesco et al.ll2008l ) and 
only the dense cores themselves were detected. However, 
the entire area will be mapped at 450 fim and 850 /im with 
SCUB A- 2 dHolland et alll200d) as p art of the Gould Belt 
survey (|Ward-Thompson et al.H2007al ). 

The 13 CO and C 18 emission in the L 1495 E and 
L 1495 SE regions mapped is very weak; in 13 CO the peak in- 
tensity towards the compact features seen in Fig.[T]measures 
only T\ 2.0 K and ^ 2.5 K, respectively. Profiles are gen- 
erally narrow (~2 km s _1 wide) and are either flat-topped 
or double-peaked (and therefore may in some regions be self- 
absorbed). 13 CO channel maps for both regions are shown in 
Figs. [2] and [3] Note the small group of compact, marginally 
red-shifted knots to the east of cores 7 and 11 in the L 1495- 
E region, and the more diffuse, slightly blue-shifted emission 
to the west and north-west. Compact features are also evi- 
dent along the south-east ridge, at both blue and red-shifted 
velocities. In both regions, however, with the possible excep- 
tion of core 14 in Fig. [3] we again see no obvious correlation 
with the H 13 CO+ cores. 

To illustrate the velocity structure observed in CO 3- 
2, we present in Figs. |4}[6] intensity- weighted radial velocity 
maps, colour-coded to show subtle changes in the centroid 
velocity of the gas (note that CO 3-2 line emission is detected 
at all locations in each image). These reveal the collimated 
blue- and red-shifted lobes of a number of outflows in L 1495, 
but also large-scale velocity gradients and "bubbles" of gas 
associated with clusters of young stars. In L 1495 E, for ex- 
ample, diffuse, red-shifted emission envelopes the cluster of 
YSOs to the east of the H 13 CO+ cores 5-8 (Fig.H the YSOs 
are marked with crosses in this figure). CO velocities to the 
west and south of this chain of cores are more blue-shifted (a 
result confirmed by the CO 1-0 observations of Goldsmith 
et al. 2008). The cluster of YSOs and the associated gas 
may be detaching itself from the rest of the cloud; alter- 
natively, the red-shifted CO may represent a bubble being 
driven eastward by the cluster of young stars. 

In L 1495 W, H 13 CO+ core 4 seems to be surrounded by 
a number of young stars; those to the north coincide with a 
region of blue-shifted CO, while those to the south coincide 
with a very distinct cloud of red-shifted gas (Fig. [5}. T Tauri 
stars are often more widely distributed than their younger 
protostellar counterparts and prestellar cores; in L 1495 E 
and L 1495 W we may therefore be witnessing the expansion 
of the natal cloud and the dispersal of these young stars. 

We also note the east-west velocity gradient seen in the 
L 1495 N region. The map in Fig. [Sjsuggests that the diffuse 
background emission to the north-west of the L 1495 bowl 
- in the unobserved regions in the centre and top of the 
figure - is likely to be red-shifted with respect to the denser 
material in the bowl. Again, this is ev ident in the CO 1-0 
observations of lGoldsmith et all l|2008l ). 

Along the SE ridge, a number of collimated outflow 
lobes (discussed further below) are superimposed onto a 
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Figure 5. Same as Fig.|4]but for the region labelled L 1495 N and L 1495 W in Fig.[T] Again, there is CO 3-2 line emission throughout 
this region. 



ridge of diffuse CO. In the centre of Fig. [6]there is some in- 
dication that the south-west edge of the ridge is red-shifted 
with respect to the north-east edge. The cores and YSOs 
in this area lie predominantly along the boundary between 
this blue- and red-shifted gas. If the cores and young stars 
delineate the high-density axis of the L 1495 ridge, then 
this distribution of blue-shifted and red-shifted gas could be 
symptomatic of rotation. 



3.2 Averaged spectra and PV diagrams 

Averaged spectra and Position- Velocity (PV) diagrams are 
shown in Figs. [7] and [8] The former illustrate the relatively 
quiescent nature of the bulk of the gas in L 1495. They also 
give an indication of the mean opacity in the three lines 
observed. If the gas is optically thin then the line inten- 
sity ratios at a given velocity should approach the abun- 
dance ratios, X[ 12 CO/ 13 CO1=70 and Xf 13 CO/ C 18 01=8.4 
(|Frerking. Langer fc Wiisonl Il982l ; I Wilson! 1 19991 ). provided 
that the excitation temperature and beam efficiency and fill- 
ing factors are the same for all three lines. Photo-dissociation 
and chemical fractionation effects in low-exti nction regions 
can also lead to enhanced isotopic ratios (jLanger et al.l 



1989). At intermediate velocities (5-8 km s _1 ), where 
line emission is detected in all three isotopologues, the 
13 CO/C ls O line ratio is indeed close to this canonical value 
(see lower panels in Fig. [7|). However, at these velocities the 
12 CO/ 13 CO ratio is of the order of 5-10, so the 12 CO emis- 
sion is clearly optically thick. For a 12 CO/ 13 CO abundance 
ratio of 70, the 12 CO opacity, T12, is given by: 

fn( 12 CO).8v _ (l- e -n 3) 

jTX( ls CO).8v (l-e-na) [ > 

where n 3 = Ti 2 /70. A 12 CO/ 13 CO line ratio of ~ 10 
thus results in T12 ~ 7.4. Again, this value will be overesti- 
mated since 12 CO is almost certainly self-ab sorbed; a lower 
abundance ratio (e.g. iGoldsmith et al] [2008) would also re- 
sult in a lower opacity. The ratio does increase towards blue- 
shifted and, particularly, red-shifted velocities. It therefore 
seems likely that, in the extended regions where the bulk of 
the outflowing gas is found, the high velocity line emission 
used to calculate the outflow parameters listed below will be 
only marginally optically thick. 

PV diagrams, some of which are shown in Fig. [S] are 
useful for distinguishing compact molecular outflows from 
turbulent motions in the large-scale cloud. The outflows are 
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Figure 7. Top: 12 CO 3-2 (blue/dot-dashed), 13 CO 3-2 
(black/dashed), and C 18 3-2 (green/full) spectra, averaged over 
the full extent of the areas mapped simultaneously in 13 CO and 
C ls O in L 1495 E and L 1495 SE. The 13 CO and C ls O spectra 
have been binned to a resolution of 0.25 km s~ 1 ; the resolution 
of the 12 CO 3-2 spectra is 0.05 km s _1 . Bottom: ratio of these 
averaged spectra; 12 CO/ 13 CO full red lines, 13 CO/C ls O dashed 
red lines. 



seen as discrete horizontal spikes; turbulence appears as re- 
gions of diffuse emission that are more extensive along the 
spatial axis (vertically in these plots) though confined to 
low velocities, i.e. to within a few km s of the ambient 
gas velocity. Indeed, we suggest that "eyeballing" outflows 
in PV diagrams is a relatively reliable way of finding them 
in HARP data cubes, since they are spatially (in one axis 
at least) as well as spectrally distinct from turbulent cloud 
motions and multiple cores seen along the same sight-line. 

We have therefore searched for outflows by stepping 
through our data cubes in latitude and then in longi- 
tude, identifying spikes in PV space that extend more than 
±2.5 km s _1 from the local ambient velocity with an in- 
tensity of >0.3 K (roughly the 2a level). A spike is defined 
as being narrower than 1' in the spatial direction (i.e. in 
longitude and/or in latitude); we adopt 1' as a conserva- 
tive upper limit to the width of an outflow lobe. Outflows 
are subsequently verified in red-shifted and blue-shifted in- 
tegrated intensity maps, most of which are presented in the 
next section. 

A number of outflows are evident in Figs. [8] In the 
L 1496 E region, for example, a compact, red-shifted knot of 
emission ~8' to the west-south- west of CoKuTau-1, which 
we later refer to as E-CO-R1, is very distinct at declination 
~28° 18' in Fig. [SJl This feature is also evident ~5' south- 
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Figure 8. CO 3-2 position-velocity diagrams measured along vertical slices through L 1495 E at RA 4:18:14.2 (a) and RA 4:17:30.1 (b), 
through L 1495 N at RA 4:16:30.2 (c), through L 1494 W at RA 4:14:22.2 (d), and through L 1495 SE at RA 4:21:11.1 (e) and RA 
4:19:28.7 (f). 



east of core 7 in Fig. [4] Similar, though less extreme spikes 
extending red-ward and blue-ward of the ambient gas are 
seen to the north of V892 Tau (not shown, though discussed 
further below). Fig. [Hp is typical of PV diagrams further 
west, toward the centre of the L 1495 bowl, where we find 
no clear-cut examples of molecular outflows, though where 
the gas generally appears more turbulent. This is particu- 
larly true in the region labelled L 1495 N in Fig. [T] where the 
emission profiles are double-peaked (top section in Fig. |HJ;). 
The CO 3-2 emission appears somewhat more ordered to the 
south and west of this region (see for example Fig. [SJ: [lower 
section] and Fig. [SJl). 

Along the south-eastern ridge in L 1495, PV spikes asso- 
ciated with high- velocity outflows are more prevalent. Exam- 
ple PV diagrams, again plotted along the declination axis, 
are presented in Fig. [8^ and[8j: these show emission from the 
complex HH 392 region and the red lobe of the collimated 
bipolar outflow associated with IRAS 04166+2706 and core 
13b, respectively (discussed further below). 



3.3 Outflows in L 1495 

When trying to identify outflows in complex regions like 
L 1495 one ideally needs a measure of the systemic velocity 
associated with each outflow driving source. This is difficult 
to do using the CO 3-2 observations alone (Fig. [SJ; even 
the 13 CO and C ls O profiles are double-peaked in some re- 
gions (Fig. [7}. We therefore use the H 13 CO + observations 
of OMK02 as a guide; these data are more likely to be op- 
tically thin and trace only the quiescent, high-density gas 
in each area. Cores are mapped throughout the bowl and 
south-east ridge in L 1495. The H 13 CO + data also suggest 
that the ambient velocity does not change drastically across 
the region (Table [TJ, a result which is supported by our CO 
3-2 observations (particularly when one scans through the 
CO 3-2 data cubes in PV space), and our more optically 
thin, though less-extensive, 13 CO and C 18 data (Figs. [2] 
and 

When constructing integrated intensity maps of the 
high- velocity blue and red-shifted gas, we therefore adopt 
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Figure 9. High- velocity CO 3-2 emission in L 1495 E, integrated 
between 4.0 and 6.0 km s" 1 (blue/thin contours) and 8.0 and 
10.0km s — 1 (red/thick contours), plotted as contours on top of 
a narrow-band H2 (+ continuum) image of the region; contour 
levels are 1.75,2,2.5,3,4,5,7,10 K km s — x . Green circles, crosses 
and triangles mark H 13 CO + cores, YSOs and HH objects, re- 
spectively. 



systemic velocities of 7.0 km s , 6.7 km s 1 and 6.6 km s 1 
for the L 1495 E, L 1495 W and L 1495 SE regions, re- 
spectively. Contour maps for select regions are presented in 
Figs. [9] - 1121 where we also over-plot the positions of the 
Taurus YSOs, known HH objects and H 13 CO+ cores for ref- 
erence. 

In these figures we make use of shallow though exten- 
sive near-infrared images of Taurus, obtained as part of 
the U.K. Infrared Deep Sky Survey currently being con- 
ducted at the United Kingdom Infrared Telescope (UKIDSS: 
iLawrence et al.ll2007l ). Images in broad-band K and narrow- 
band H2 2.122 fim emission (hereafter referred to simply as 
H2) have been secured for the entire Taurus- Auriga-Perseus 
complex. The wide-field camera (WFCAM) used to obtain 
these data, the data reduction proced ure, and th e WF - 
CA M archive ar e desc ribed in detail bv lLucas et al.l (|2008l ) 
and lDavis et alj (|2008l ). Note that we also use the acronym 
"MHO" for the molecular hydrogen emission-line objects - 
the shock-excited infrared counterparts to HH objects - ob- 
served in L 149E0. A complete cata logue of all known g alactic 
MHOs is available on-line0 and in lDavis et~afl (^Old ). 

Below we discuss the star forming regions in L 1495 
individually. 



3.3.1 L 1495 E: CoKu Tau-1 and HH 156 

L 1495 E harbours a cluster of about a dozen TTSs, many 
of which have been identified throu gh X-ray observation s 
(the V 410 X-ray and Anon# sources; Istrom fc StroirJl994h . 
OMK02 label five cores in this region. In the integrated in- 
tensity maps in Fig. [T] the CO is diffuse and optically thick 
across much of the region; only a few compact peaks are 
detected, all of which are marginally red-shifted and unre- 
lated to the H 13 CO + cores (Fig. though note that the 
H 13 CO + map does not cover the region south of the young 
star Xray 7 and east of core 8). 

Two seemingly unrelated red-shifted knots are identi- 
fied as molecular outflow features, E-CO-R1 and E-CO-R2. 
Neither is obviously related to any of the YSOs or H 13 CO + 
cores; their progenitors remain unknown, although candi- 
date sources surround E-CO-R1, and this object does coin- 
cide with a peak in our 13 CO integrated intensity map in 
Fig. [T] Note also that in the H2 images we detect no line- 
emission features in this region. 

In Fig. [9] CoKu Tau-1 seems to be associated with a 
third compact red-shifted CO feature. However, this low- 
velocity knot of emission fails our search criteria for out- 
flows (described above) and so is not identified as such here. 
CoKu Tau-1 does power HH 156 A and B, two compact HH 
objects offset ~2" and ~12" to the south-south-west; the 
CO peak could represent the counter-lobe, though deeper 
and/or higher-resolution CO data are needed to be sure. 
(The lack of blue-shifted CO around HH 156 may be a re- 
sult of the HH flow exiting the cloud, or may be due to the 
dispersal of the core around this TTS.) 

Apart from CoKu Tau-1, none of the young stars in the 
region are obviously driving CO outflows; presumably most 
are too evolved. Of the dozen YSOs to the east of cores 5-8 in 
Fig- El all bar three have neutral Spitzer-TRAC colours ([3.6]- 
[4.5]~0; [5.8]-[8.0]~0), consistent with relatively evolved 
(weak-line) TTSs and an absence of mid-infrared exces s due 
to a lack of circumstellar material |Luhman et "ai]|2006l . note 
that V 892 Tau is saturated). Only CoKu Tau-1, CZ TauA 
and DD Tau A have 2MASS near-IR and Spitzer mid-IR 
magnitudes that brighten towards longer wavelengths, con- 
sistent with them being protostars or very young T Tauri 
stars and therefore possible molecular outflow candidates. 
This low number of protostars in L 1495 E is consistent 
with the paucity of outflows in our H2 and CO data. 



3.3.2 L 1495 W: CW Tau and HH 220/826-828 

The high-velocity CO in L 1495 W is shown in Fig. ITP1 
Again, the positions of candidate YSOs a re marked, al- 
thoug h note that the Spitzer observations of iLuhman et all 
(2006) do not cover this region. The most dramatic feature 
in this area is the bubble of red-shifted gas that surrounds 
a small group of young stars, YSOs l,2,3|j and CIDA-1 (see 
also Fig. [5} . More compact blue-shifted features are also ob- 
served in L 1495 W. 

At least four molecular outflows exist in this region: 



2 The acronym MHO has recently been approved by the Interna- 
tional Astronomical Union (IAU) Working Group on Designations 
and has been entered into the Dictionary of Nomenclature of Ce- 
lestial Objects (http://cdsweb.u-strasbg.fr/cgi-bin/Dic7MHO) 

3 http://www.jach.hawaii.edu/UKIRT/MHCat/ 



4 Kenyon et al. (2009) refer to these objects as MHO 1, MHO 2 
and MHO 3, using the acronym MHO for "Mill House Observa- 
tory". We identify these sources with YSO so as to distinguish 
them from Molecular Hydrogen emission line Objects (MHOs). 
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Figure 10. Same as Fig. [9]but for the L 1495 W region. Left - high-velocity CO 3-2 integrated between 3.7 and 5.7km s — 1 (blue/thin 
contours) and 7.7 and 9.7km s _1 (red/thick contours). Contour levels are 1.0 (blue-shifted only), 1.2, 1.5, 2, 3, 4, 5, 7, 10 K km s _1 . Right - a 
close-up view of the region to the south-east of core 4 (marked with a dashed square at left); contours measure 1,2,3,4,5,7,10 K km s _1 . 



(1) a compact north-south bipolar flow, W-CO-flowl, cen- 
tred ~1' south of (and seemingly unrelated to) H 13 CO + 
core 4, (2) a collimated blue-shifted lobe extending north- 
westward from YSOs 1 and 2 (which presumably also drive 
some [or all] of the red-shifted CO to the south-east), (3) 
a compact blue-shifted CO feature, W-CO-B1, associated 
with HH 828, and (4) a second blue CO knot, W-CO-B2, 
to the south-west of core 4, which may well be driven by 
IRAS 04108+2803A; this CO knot is certainly extended to- 
wards this object, which is also rather red (J-H = 3.10, H-K 3 
= 2.32). 

Two additional flows are detected i n HH and H2 emis- 
sion HH 220/826, a jet from CW Tau jMcGroartv fc Ravi 
l2004h , and an arc of HH/H 2 emission, HH 827/MHO 715, 
which may be associated with low-velocity red-shifted CO 
extending southward from CIDA-1 (though again this emis- 
sion fails our outflow criteria). 

YSO 1/2, YSO 3 and CIDA-1 have near-IR 2MASS 
colours consistent with youth. The TTSs CW Tau, 
V773 Tau A and FM Tau are more evolved and appear to be 
associated with a completely separate region of star forma- 
tion. CW Tau is thought to be the most active classical TTS 
in this region. Even so, we detect no H2 line emission from 
its HH objects, and only HH 828 appears to be associated 
with high-velocity (blue-shifted) CO emission. 



3.3.3 L 1495 SE: HH 390-391 

The lower portion of the L 1495 ridge is populated with 
a number of young stars and HH flows. The area around 
cores 13a, 13b, 14 and 15 is shown in Fig. 1111 Two of these 
four H 13 CO + cores (13b and 14) are associated with IRAS 
sources and HH objects. 

IRAS 04169+2702 (core 14) drives a collimated, pos- 
sibly precessing bipolar molecular outflow and HH object 
(HH 391 A); the H H object is also iden tified in H2 emission, 
MHO 731 (see also lGomez et al.lll997l ). 



A second, rather spectacular bipolar CO outflow em- 
anates from core 13b. This flow is particularly striking 
in Fig. H The source of the outflo w, IRAS 04166+2706, 
is ide ntified at submm wavelengths (|Santiago-Garcfa et al.l 
2009), though it does not appear in the 2MASS point 
source cat alogue nor in the Y SO lists of iLuhman et al.l 
( 2006) and iKenvon et all (|2009l V It must therefore be par- 
ticularly young, as expected for molecular outflow and 
MHO progenitors (Hatchell et al. 2007b; Davis et al. 
2008, 2009). Very faint H 2 features, MHO 717 A/B and 
MHO 718, are observed in both flow lobes. This outflow 
is wel l-known from previous CO observations: iTafalla et al.l 
|2004h report the discovery of a collimated, Extremely High 
Velocity (EHV) outflow fr om IRAS 04166+2706, which 
ISantiago-Garci'a et all (|20Qgh have since observed in SiO 2-1 
emission. Both groups present CO 2-1 spectra with striking 
peaks in the line wings at -35 km s _1 and +45 km s _1 . 
However, this EHV outflow component is only marginally 
detected (at the 2-3cr level) in our 3-2 data. Based on a 
noise level of TX(rms) ~ 0.1 K in our data (equivalent to 
a main beam brightness temperature of ~ 0.2 K), we set 
an upper limit of ~0.5 for the 3-2/2-1 intensity ratio (note 
that the beam sizes are comparable). This ratio suggests 
that the clumps or "molecular bullets" associated with the 
EHV flow from IRAS 04166+2706 ar e relatively cold and/or 
diffuse (n H2 < 10 4 cm" 3 , T ^ 50 K; lYeh et al]|200Sl ). 

Core 13a drives a much less well defined molecular 
outflow. MHO 716 may be associated with a red-shifted 
outflo w lobe; HH 390 (w hich extends towards the south- 
west; I Gomez et al l |1997T > would then be excited in the 
blue lobe. The 2MA SS s ource J04194148+27 16070, which 
ILuhman et"afl (|2006h and lKenvon et al.l (|20Qgh mis-identify 
as IRAS 04166+2706, is situated ~50" to the north-east of 
core 13a, and is a possible driving source for this outflow. 

Considering these three cores (13a, 13b and 14) and 
their associated YSOs together, IRAS 04166+2706 is pre- 
sumably the youngest of the three; this is consistent with 
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Figure 11. Same as Fig.[9]but for the region around cores 13-15 in L 1495 SE. Top-left - high-velocity CO 3-2 integrated between 3.6 
and 5.6km s _1 (blue/thin contours) and 7.6 and 9.6km s _1 (red/thick contours); levels are 2,3,4,5,7,10 K km s — 1 for both blue- and 
red-shifted emission. Right and Bottom-left - close-up views of the outflows associated with cores 13a, 13b and 14: contour levels are 
1.5(red only), 2, 3, 4, 5, 7,10 K km s _1 in the former and 2,3,4,7,10 K km s _1 in the latter. 




Figure 12. Same as Fig. |9]but for the region around cores 16a, 16b and 17 in L 1495 SE. Left - high-velocity CO 3-2 integrated between 
3.6 and 5.6km s — 1 (blue/thin contours) and 7.6 and 9.6km s _1 (red/thick contours); levels are 1,2,3,4,5,7,10 K km s" 1 for both blue- 
and red-shifted emission. Right - a close-up of the outflows associated with cores 16a and 16b; contour levels are 1,2,4,7,10 K km s . 
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its powerful EHV CO outflow, its H2 shocks, and the ab- 
sence of HH objects. IRAS 04169+2702 is probably some- 
what more evolved, having infrared colours that are consis- 
tent with a Class I protostar (H-K a = 2.6; [3.6]-[5.8] = 2.1). 
J04194148+2716070 is bluer and is almost certainly a TTS, 
which would explain its association with an HH object and 
the lack of well-defined CO 3-2 outflow lobes. 

Further south, HH 391 B and C represent a chain of 
optical f eatures that extend over ~1' in a north-south di- 
rection |Gomez et al.l fl997l ). These objects coincide with 
large, complex bubbles of both blue- and red-shifted emis- 
sion. HH 391 B/C are unlikely to be associated with the 
IRAS 04169+2702/HH 391 A outflow. No H 2 emission was 
detected in this region, and the relationship between these 
shock features and the complex CO emission in this area is 
unclear, although high- velocity CO does coincide with these 
objects (we label the more clearly-defined red-shifted emis- 
sion SE-CO-R1). Moreover, a 13 CO peak does coincide with 
HH 391 B/C in our integrated intensity map in Fig. [TJ which 
may be a sign of youth. 

8.3.4 L H95 SE: FS Tau A/B (Haro 6-5B) and 
HH 157/276/392 

Moving further south-east along the L 1495 ridge, molecu- 
lar outflows are detected around IRAS 04181+2654 and the 
well-studied FS Tau A/B (Haro 6-5B) system (Fig. [12]). 

IRAS 04181+2654 A is closest to core 16b and ap- 
pears to drive a bipolar CO outflow in a roughly east- 
west direction. The blue lobe of this flow may overlap 
with a second bipolar outflow, SE-CO-flowl, which drives 
the compact shock features HH 392 C/MHO 719. There 
is no source candidate for SE-C O-flowl in the 2MAS S or 
IRAS point source ca talogues, and lLuhman et all (|2006T l and 
iKenvon et all ||2009h do not identify a YSO near these ob- 
jects. IRAS 04181+2654 A is resolved by 2MASS and ex- 
hibits an excess consistent with a protostar (J-H = 3.58, 
H-K s = 2.30). 

IRAS 04181+2654 B lies ~30" north-west of star A; 
this nebulous infrared source is associated with an arc of 
HH emission, HH 392 B (tentatively detected here in H2) 
and also possesses considerable excess in the near-IR (J- 
H = 5.11, H-K s = 2.70). The third source in the region, 
2MASS source J042 10795+2702204, has less extreme near- 
IR colours (J-H = 1.8, H-K s = 1.5), though may nonetheless 
also contribute to the complex high-velocity CO evident in 
Fig. 1121 The "kidney" shape of the blue- and red-shifted 
features coincident with HH 392 A and B certainly suggest 
the presence of multiple flow lobes. 

Approximately 1.5' to the north-east of 
IRAS 04181+2654 there is a faint knot of H 2 emis- 
sion, MHO 720, which coincides with a knot of red-shifted 
CO, SE-CO-R2, and a very marginal detection of blue- 
shifted CO offset ~15" to the south-west. In all, there are 
at least four molecular outflows within a 5' radius of cores 
16a and 16b. 

FS Tau A (J-H = 1.47, H-K s = 1.73) and FS Tau B (J-H 
= 1.06, H-K s = 1.60; also known as Haro 6-5B) are a pair of 
TTSs, the latter being associ ated with the spectacula r HH 
jet and bow shock HH 157 (jEisloffel fc Mundtl Il998f l. We 
detect red-shifted CO emission from this outflow, though no 
H2 line emission. HH 276 is a chain of faint HH knots that 



crosses the HH 157 jet ~1' to the north-east of FS Tau B, 
although no high-velocity CO is clearly identified from this 
flow. 

Lastly, south-east of FS Tau A, we note the very tenta- 
tive detection of a bipolar CO outflow from CFHT-21, orien- 
tated NE-SW, roughly parallel with the FS Tau B jet. This 
marginal detection is consistent with the near-IR colours of 
this TTS (J-H = 1.5, H-K s = 1.0). SE-CO-B1 is a com- 
pact though massive high-velocity knot with a tail extend- 
ing along the south-east ridge. This feature is, however, not 
related to any known YSOs in the region. 



3.4 Outflow parameters 

The physical parameters for the molecular outflows de- 
scribed above are listed in Table [2] (we also list HH out- 
flows not detected in high- velocity CO emission). The mass 
in the molecular outflow lobes (the mean of both lobes if 
the outflow is bipolar) is derived from the column density 
integrated across the extent of the high- velocity line wings, 
which in turn is calculated from the integrated antenna tem- 
perature, corrected for the main-beam efficiency (we assume 
a beam filling factor of unity), and an estimate of the CO 
3-2 excitation temperature. 

The excitation temperature, T ex , in the ambient gas can 
be estimated from the peak temperature in the CO 3-2 line 
profiles, pr ovided the gas is optically thick though not self- 
absor bed (|Pineda. Caselli fc Goodman! 120081 ; iBuckle et al.l 
120101 ). However, T ex in the ambient gas is not necessarily 
the same as T ex in the entrained outflow gas, where the ki- 
netic temperature may be higher and the gas is almost cer- 
tainly compressed as it is s wept up by the under l ying jet and 
HH/MHO bow shocks (e.g. IHatchell et al.ll 19991; Davis et all 



2000al; ICiannini. Nisini fc Lorenzettill200ll ; iLee et al.ll2002l : 



Kempen et al.ll2009h . We therefore adopt a somewhat 



higher value than is implied by (a) the average CO 3- 
2 spectra in Fig. (TJ ~ 3 K; T ex ~ UK), which 
include the more diffuse regions away from the central 
south-east ridge and star forming cores, or (b) the bright- 
est spectra observed in L 1495 (T4 ~ 17 K; T ex ~ 
35 Kn. We instead use a value of 50 K, consistent with 
the range in excitation temperature derived from CO 
line ratios in oth er outflows - albeit assuming optical! 
thin emission (e.g. IDavis. Smith fc Moriartv-Schievenlll9"9T 



IHatchell et~al]|l998l : IDavis et alj|2000al: Ivan Kempen et al.l 
2009). The mass is in any case relatively insensitive to tem- 
perature, varying by only 40% o ver the temperature range 
20-100 K (IHatchell et al. 2007a). Assuming an abundance 
of A cq = 10~ 4 (|Frerking. Langer fc Wilson! 1 19821 ; IWilsonl 
1999), the column density (in units of cm -2 ) is then given by 
Nh2 = 2.5 x 10 19 ( T£dv), where the integrated line wings 
are in units of K km s" 1 (|Hatchell et al.ll2007al ). J T£dv is 
measured from the integrated high- velocity blue-shifted and 
red-shifted outflow maps used to plot the contours in Figs. [5]- 
1121 the integrated flux measured in an ellipse that envelopes 



5 In an isothermal slab, assuming local thermodynamic equilib- 
rium and optically thick emission, the excitation temperature, 
T ex , is related to T(,, the line peak main beam brightness temper- 
ature, by: T e3; (3 - 2) = 16.59/Zn[l + 16.59/(T b ( 12 C*0) + 0.036)] 
llPineda. Caselli fc Goodmanll2008l) . 
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Table 2. Parameters for the outflows in L 1495. The top nine flows are in the bowl; the remaining flows are in the south-east ridge. 



YSO or Outflow* 1 


Arca a 


YSO 
type a 


RA a 
(J2000.0) 


Dec a 
(J2000.0) 


CO flow b 


HH/MHO c H 13 CO + 


M d 
(Mq) 


km s- 1 


(Mq 

km 8-1) 


E d 
(X10 34 
J) 


CoKu Tau-1 


E 


TTS 


4:18:51.5 


28:20:26 


No? 


HH156 










E-CO-Rl e 


E 




4:18:39.5 


20:20:20 


Red 




0.0098 


2.3 


0.022 


5 


E-CO-R2 e 


E 




4:18:15.0 


28:18:30 


Red 




0.0256 


4.0 


0.102 


41 



YSO 1/2 

CIDA-1 

CW Tau 

W-CO-Bl e 

W-CO-flowl e 

IRAS 04108+2803A? 



Class I 
TTS 
TTS 



4:14:26.3 
4:14:17.6 
4:14:17.0 
4:14:10.7 
4:14:12.0 
4:13:57.4 



1:06:03 
1:06:10 
1:10:58 
1:14:40 
1:08:30 
1:09:10 



HH827/MH0715 

HH220/826 

HH828 



0.0028 
0.0038 
0.0055 



3.1 
3.1 
3.3 



0.009 
0.012 
0.018 



J 041 941 484- 2 71 60 70 


SE 


TTS 


4:19:41.5 


27:16:07 


No 


HH390/MHO716 


13a? 










IRAS 04166+2706-f 


SE 


Class 0? 


4:19:42.4 


27:13:24 


Bipolar 


MH0717/71S 


13b 


0.0217 


8.9 


0.193 


172 


IRAS 041694-2702 


SE 


Class I 


4:19:58.5 


27:09:57 


Red 


HH391A/MH0731 


14 


0.0301 


6.1 


0.184 


112 


SE-CO-R1 


SE 




4:19:52.5 


27:05:30 


Red? 


HH391B/C 




0.0592 


6.2 


0.367 


228 


SE-CO-R2 


SE 




4:21:17.1 


27:02:50 


Red 


MHO720 




0.0021 


2.6 


0.005 


1 


J042 10795 + 2702204 


SE 


Early TTS 


4:21:08.0 


27:02:20 


Bipolar 


HH392A 




0.0157 


5.4 


0.085 


46 


IRAS 04181 + 2654A 


SE 


Class I 


4:21:11.5 


27:01:09 


Red? 




16b 


0.0049 


2.6 


0.013 


3 


IRAS 04181 + 2654B 


SE 


Class I 


4:21:10.4 


27:01:37 


? 


HH392B 


16b? 










SE-CO-flowl e 


SE 




4:20:56.0 


27:00:00 


Bipolar 


HH392C/MH0719 




0.0153 


3.2 


0.049 


16 


FS Tau B 


SE 


Early TTS 


4:22:00.7 


26:57:32 


Red 


HH157 




0.0209 


4.9 


0.102 


50 


CFHT-21 9 


SE 


TTS 


4:22:16.8 


26:54:57 


Bipolar? 














SE-CO-Bl e 


SE 




4:22:22.0 


26:48:00 


Blue 






0.0251 


3.8 


0.095 


36 



a Thc most likely driving source of the CO outflow, the area in which it is found in L 1495, its YSO type (based on near- and/or mid-IR colours) and its coordinates. 
If no YSO source is identified the label used for the CO outflow or high- velocity CO peak is given. The coordinates then refer to the approximate midpoint between 
the bipolar lobes, or the peak in the CO emission if the flow is monopolar. 

b Indication of whether the CO outflow is bipolar or monopolar (with a red-shifted or blue-shifted lobe clearly defined). "No" in this column indicates that no CO 
outflow was detected. 

c Associated HH objects, MHOs and H 13 CO+ core (from OMK02). 

d Thc mass, maximum radial velocity (the difference between the nominal ambient velocity and the velocity in the line wings at the 2ct noise level); the momentum 
and the kinetic energy of each outflow lobe. Mean values arc quoted if the flow is bipolar 

CO outflows with no clearly identified driving source. 
f Driving source not listed in the YSO catalogue of |Luhman et alj <2u06j ), [^enyon et "al~| J20091 ), or detected by 2MASS. 
^CO outflow only a marginal detection; requires confirmation with deeper observations. 



each outflow lobe is "sky subtracted" using the flux mea- 
sured in an outer annulus (this eliminates diffuse red-shifted 
or blue-shifted gas not associated with the outflow); the de- 
rived column density is then scaled by the HARP beam size 
and the factor 1.3xrriH2 (where mm is the mass of an H2 
molecule) to give the mass in each outflow lobe. 



We assume that the wings of the CO 3-2 emission are 
optically thin, so the values in Table [2] will be lower limits: 
the mass, M, momentum, P, and kinetic e nergy, E, are prob- 
ably underestimated by a f actor of 2-4 (jCabrit fc Bertoutl 
1990; lHatchell et all 12007 j ). Furthermore, the radial ve- 
locity, \v — v\ (v is the velocity in the line wings where 
the emission reaches the 2a noise level; v is the ambient 
velocity from Table [T]), used to calculate the momentum 
and energy, will also be underestimated, by a factor 1/sinz, 
where i is the inclination angle with respect to the plane 
of the sky. Obviously for flows close to the plane of the 
sky this factor can be quite considerable, although flows 
with i 3> 60° will be difficult to detect in our CO ob- 
servations; for an inclination angle of 60° P and E will 
be underestimated by an additional factor of 2 and 4, re- 
spectively. Even so, the parameters listed in Table [2] are 
not unusual for outflows from low-mass young stars. They 
are marginally lower than (within an order of magnitude 
of) those typically observed for low-mass YSO outfl ows in 



Orion, Ophiuchus and Perseus (e.g . [Knee & S andelll 



Williams. Plambeck fc Heverl 120031; lArce fc Sargentl 



2000; 



2006; 



0.0 
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Figure 13. Near-IR colour-colour diagram for the HH jet and 
molecular outflow sources in L 1495, which are marked with blue 
squares and green diamonds, respectively. The black dots lower- 
left represent the l ocus of intrinsic colours for main sequence 
dwarfs and giants jKoornnedll983h : the dashed line indicates 
the locus of T Tauri stars (Meyer, Calvet & Hillenbrand 1997). 
The full parallel lines arc redden ing vectors for the main sequence 
stars up to A v = 30 mag iRieke fe Lebofskvll 19851) . 



Bussmann et al.ll2007l ; IStanke fc Williamdl2007f l. 
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4 DISCUSSION 

4.1 Using outflows to distinguish starless from 
protostellar cores in L 1495 

L 1495 contains the largest concentration of young stars in 
the Taurus region. The majority of these are found in the 
L 1495 bowl, particularly the denser eastern half of the bowl 
(the molecular gas in the western bowl region is more ex - 
tensive and generally more diffuse; G oldsmith et al. 20081) 
and in the lower half of t he south-east ridge ( Luhman et al.l 
l2006l ; lKenvon et al.llioodl '). The low number-density of YSOs 
in the upper, north-western half of the ridge, where OMK02 
nevertheless find a number of massive, dense cores (9a, 9b, 
10a, 10b and 12), suggests that star formation has yet to 
take place in this area. The absence of HH objects, MHOs, 
and CO outflows in this region supports this interpretation. 

The molecular outflows in L 1495 are listed in Tabled 
along with associated HH objects, MHOs, H 13 CO + cores 
and candidate YSO driving sources. YSO source classifica- 
tions are based on near-IR and/or mid-IR photometry (as 
described in section 1). In Fig. [13] we plot the 2MASS colours 
of the HH jet and outflow sources in Table Hon a near-IR 
colour-colour (CC) diagram. As expected the majority of the 
sources lie to the right of the reddening band (the two paral- 
lel lines), although notably the sources with no CO outflow 
all lie in the region associated with TTSs, while the young 
stars that do drive molecular outflows are all much redder (in 
both J-H and H-K a ). Clearly the molecular outflow-driving 
sources are more embedded than their HH jet-driving coun- 
terparts. 

Of the 22 cores labelled in Fig. [T] only four - 13a, 13b, 
14 and 16b - appear to be associated with YSOs (in each 
case the YSO is located within the 50% integrated intensity 
contour in the H 13 CO + maps). This relatively small frac- 
tion of cores with YSOs is consistent with the rest of Tau- 
rus; OMK02 find that only 22% of their condensations are 
associated with embedded sources. In L 1495 these "cores 
harbouring stars" certainly seem to be protostellar; at least 
three of the four cores (13b, 14 and 16b) are associated with 
molecular outflows; the CO flow from 13a is a marginal de- 
tection, though this core is also probably associated with 
an HH jet a nd MHO. As has been noted in other outflow 
surveys (e.g. iDavis et al.ll2008l ; lHatchell fc Dunham! 120091 ) , 
this suggests that CO outflows are useful tracers of the lo- 
cations of protostellar cores. Indeed, outflow surveys may 
be used in combination with near- and mid-IR photometric 
studies of star forming regions to establish more accurately 
the fraction of cores that do harbour accreting protostars. 

Spitzer, when combined with near-IR data from e.g. 
2MASS or WFCAM, can be a powerful tool for searching 
for Class 0/1 protostars (although saturation can be a prob- 
lem with Spitzer data, particularly at longer wavelengths). 
This is especially true in nearby, low mass star forming re- 
gions like_Tairrus, where extinction and crowding are mini- 
mal l|Evans et alj|2009h . A few recent studies have revealed 
the presence of extremely faint protostellar objects inside a 
handful of cores that were previously t hought to be starless 
l|Crapsi et al.ll2005l ; iBourke et al.ll2006l ). But such low-mass 
cores may be below the resolution and sensitivity limits of 
theH 13 CO+ survey discussed here, and these objects are un- 
likely to drive powerful CO outflows or bright H2 jets. The 
YSO list, derived from optical, near-IR and mid-IR pho- 



tometry and spectroscopy (jKenvon et alJ koOQ). is therefore 
expected to be relatively complete for the core mass range 
under scrutiny here. 

A number of high-velocity CO features appear to be 
without driving sources in Tabled In some cases their pro- 
genitors may be nearby; however, in a few (W-CO-flowl, 
SE-CO-flowl, and SE-CO-R1/HH 391), the location of the 
driving source is clearly defined, yet they are still unde- 
tected. These flows are probably driven by Class sources 
which require deeper mid-IR or far-IR/submm observations 
i|j0rgensen et al.ll2007l.l200Sl ). 

From our comparison of the H 13 CO + cores mapped by 
OMK02 to the YSO source list and outflows mapped here, 
at first sight the ratio of starless cores to protostellar cores 
in L 1495 seems to be high; as noted earlier, of the 22 cores 
located in L 1495, only 4 appear to be coincident with known 
YSOs (although all are associated with outflows). 

However, two cores, 3 and 4 , are found in a region that 
was not observed with Spitzer. iLuhman et alj |2006) esti- 
mate that pre-Spitzer studies of the YSO population in Tau- 
rus are only complete to about 80% . They also point out 
that Spitzer is very well suited to finding disk-bearing stars 
down to masses of the order of 0.01 Mq through consid- 
erable amounts of extinction (A v ~ 100). The observations 
of OMK02 are sensitive to core masses in the range 3.5 Mq 
< M < 20.1 Mq. Even for a very low star formation effi- 
cienc y (low protostar-to-core mass ratio) of just a few per- 
cent (jWard-Thompson et al ]|2007bh . if these "starless" cores 
did contain protostars Spitzer would detect them. Given the 
small number statistics here, it is therefore possible that 
sensitive mid-IR observations could uncover protostars as- 
sociated with both cores 3 and 4, in which case six of the 22 
cores would be associated with protostars. (Note also that 
core 4 may be associated with the CO outflow W-CO-flowl.) 

Hence, in L 1495 it seems reasonable to assume that we 
have 5±1 protostellar cores and 17±1 starless or prestellar 
cores; the ratio of cores that are associated with young stars 
to those that are not is in the range 0.2-0.4. Bu t how many of 
the starless cores are prestellar? OMK02 and[M 

izuno et al.l 

(1994) before them suggest that the fraction of H 13 CO + 
cores in L 1495 that exceed the Jeans mass, are gravita- 
tionally bound, are undergoing some form of collapse, and 
therefore are prestellar in nature, is probably high: note that 
theH 13 CO+ cores identified by OMK02 are compact (radius 
^ 0.1 pc) and of high density (^ 10 5 cm -3 , the critical den- 
sity for excitation of the H 13 CO + line used in their study). 
They find that 13 of the 22 cores (60%) in L 1495 have an 
H 13 CO + column density mass that exceeds the virial mass 
(these 13 include cores 3, 4, 13a, 13b and 16, i.e. five of 
the six candidate prestellar cores). This leads us to believe 
that 9 of the 22 are neither prestellar nor protostellar, so 
we identify these as being starless. The ratio of prestellar to 
protostellar cores in L 1495 is therefore in the range ~1.3- 
2.3, while the ratio of starless to prestellar cores is in the 
range ~1— 1.3. 

Although based on modest statistics, the above ar- 
guments indicate that the prestellar phase is about as 
long lived as the protostellar phase, and that in L 1495 
there are as many prestellar cores as starle ss cores. Sim- 
ilar re s ults have been f ound in Perseus (Hatchcll et al.l 
2007al; |j0r~t rensen et al. I 120071 ; lHatchell fc Fulled 120081 : 



Davis et al.ll2008l ) and in other low mass star forming re- 
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Figure 14. (a) Plot of H 13 CO + core mass against core radius; (b) plot of H 13 CO + line width against line intensity measured towards 
the peak in each core; (c) and (d) H 13 CO" 1 " core molecular gas density and core radius plotted against the mass of the associated outflow. 
All core parameters are taken from Table 2 of OMK02. In each panel, cores are marked with red dots and are labelled according to 
OMK02; cores associated with YSOs are marked with green diamonds; cores associated with molecular outflows or HH jets are marked 
with open blue squares (note that we were unable to measure the molecular flow mass for core 13a). 



gions l|Visser. Richer fc Chandle rl [2003 ; lEnoch et al.1 2007; 
lEvans et al.ll2009h . 

lEvans et al] {2009) find the Class lifetime to be 
~10 5 yrs and the Class and Class I lifetimes combined 
to be about 5 x 10 5 yrs. It is worth mentioning that these 
timescales are more-or-less consistent with the range in dy- 
namical ages, 10 4 - 10 5 yrs, measured for molecular outflows 
(Arce et al. 2007, and references therein); this is to be ex- 
pected, given the statistical evidence that most protostars, 
partic ularly those associ ated with cores, drive molecular out- 
flows l|Davis et al.ll2009l ). 

lEvans et al.l (|2009T ) also find the prestellar core lifetime 
to be of the order of 5xl0 5 yr. OMK02, from their analysis 
of the data discussed here, estimate a very similar timescale 
(~ 4xl0 5 yrs) for starless condensations in Taurus. If we 
assume a value of 5xl0 5 yrs for the protostellar lifetime 
of our sources, then we arrive at a prestellar lifetime of 
^9(±3) x 10 yrs. This is slightly longer than the timescale 
found by Evans et al., although it is still consistent to within 
the combined errors. Furthermore, we find the lifetime of 
starless cores that are not gravitationally bound to also be 
of the order of 10 6 yrs. 



4.2 The relationship between cores and outflows 

OMK02 estimate various parameters for the H 13 CO + con- 
densations in Taurus, including size, density and mass. We 
plot a number of these in Fig. 1141 in Fig. I14f a~) we chart core 
mass based on the integrated intensity in H 13 CO + emission 
against core radius (see OMK02 for details); a plot using 
the core virial mass, also calculated by OMK02, looks very 
similar to this figure. In Fig. Ob) we plot the H 13 CO + 
line width against H 13 CO + peak antenna temperature; both 
have been measured by OMK02 at the peak position in each 
of their H 13 CO+ cores. Cores associated with embedded 
stars and outflows or HH jets are identified in both figures 
with green diamonds and open blue squares, respectively. 

The plots in Fig. I14f a) and (b) suggest that the pro- 
tostellar, outflow-driving cores are generally less massive, 
more compact and more quiescent than many of the ob- 
served condensations (although clearly not all of the com- 
pact H 13 CO + cores are associated with outflows). At first 
glance this seems to run against expectation, since other 
studies have found that higher mass core s are more likely to 
be p restellar (e.g. |j0rgensen et al.1 120071 ; lHatchell fc Fuller! 
2008). However, one should note that the more massive 
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H 13 CO + cores are often larger (Fig. I14f a)) and less dense 
than their less massive counterparts; these objects may in 
fact harbour multiple, unresolved cores. Alternatively, these 
massive H 13 CO + condensations may be younger than the 
protostellar cores, or may simply be gravitationally unbound 
and therefore not prestellar in nature. 

For the H 13 CO+ cores that are associated with out- 
flows and protostars, we find no correlation between core 
mass and outflow mass or kinetic energy (plots not shown). 
If anything, the more massive cores (4 and I6b) seem to 
be associated with the least massive flows. However, this 
comparison does not take into account the size of each core; 
Fig- I14f c) and (d) indicate a possible correlation between 
core density and outflow mass, and an inverse correlation 
between core size and outflow mass, albeit for a very small 
number of cores and outflows. The implication here is that 
the build-up of outflow mass increases as the protostellar 
core contracts. This idea clearly requires further investiga- 
tion, both observationally and theoretically. Even so, our 
results are consistent with th e early analysis of the H 13 CO + 
data by iMizuno et all l| 19941 ). who find that condensations 
with embedded stars in Taurus tend to be smaller in size, 
more dense, and less extended, than those without. Simi- 
lar results were found i n the submm continuum survey of 
IWard- Thompson et all l|l994h . 

Finally, we have c onducted a statis tical analysis of the 
data plotted in Fig. ll4l (I sobe et al.1 1990). This suggests that 
there is a 93% chance that the correlation between core mass 
and core radius is significant (i.e. there is a 7% chance that 
the sample is random), for cores with stars and for cores 
with outflows; for the entire sample there is a much higher 
chance (>99%) that the correlation is real, in part because 
of the much larger number of data points. Likewise, we find 
a >99% chance that for all cores the line intensity is related 
to the line width (Fig. 114( b)). although when only cores with 
stars or cores with outflows are considered the distribution 
of points is far less likely to be real (—65% and <50%, re- 
spectively). There is a 98.5% chance that the linear relation- 
ship between outflow mass and core density in Figs, 114( c) 
is real, though a reduced likelyhood (85%) that the distri- 
bution between outflow mass and core mass in Figs. 114( d) 
is significant. Strictly speaking, only the linear distribution 
of points in Figs. 114( c) is potentially real, although clearly 
both plots would benefit from additional data. 



4.3 Star Formation Efficiency in L 1495 based on 
outflow statistics 

In their study of outflow activity in Perseus- West, 
I Davis et alj (|2008l ) note that the more massive clouds are as- 
sociated with a greater number of outflows. They use cloud 
masses derived from extinction mapping by Kirk, Johnstone 
& Di Francesco (2006), who note that this is arguably the 
best way of tracing the diffuse atomic and molecular gas 
(A v ~ 3 — 7) in large molecular clouds. Kirk et al. also 
stress that only —5% of the cloud mass resides in regions 
with H 2 column densities > 5 x 10 21 cm -3 (A v > 2.5): 
in other words, only 5% of the gas mass is locked in the 
starless and protostellar cores mapped in H 13 CO + and in 
submm dust continuum emission with, for example, SCUBA 
or SCUBA-2. Based on these measurements, in Perseus- 
West Davis et al. estimate that there is one outflow for 



every 44-88 Mq of ambient material. This range is not in- 
consistent with a canonical value for the star formation ef- 
ficiency (SFE = My so/\M Y so + M corc ]) of -10-15% (e.g. 
Ijprgensen et al.ll2007T ). if the cloud-to-core mass ratio is in- 
deed —20, and the protostars have masses of the order of 

~0.5 Mq. 

From their large-scale CO 2-1 an d 13 CO 2-1 observa- 
tions of Taurus. [Goldsmith et ail (|2008h find that the L 1495 
bowl and south-east ridge (regions they label L 1495 and 
B 213) have masses of the order of 2.6 x 10 3 Mq and 
1.1 x 10 3 Mq, respectively. Their observations are sensi- 
tive to H2 column densities in excess of — 1 x 10 21 cm -3 
(A v > 0.5), so somewhat similar to the columns of molecu- 
lar gas probed by the extinction mapping technique used 
in Perseus. In L 1495 we list (in Table 0) nine outflows 
in the bowl and 12 in the L 1495 ridge; we therefore es- 
timate roughly one outflow per 289 Mq in the bowl, and 
one per 92 Mq of gas in the ridge. These values are some- 
what higher than in Perseus, especially in the bowl region. 
This is likely due to the more evolved state of this region, 
where the number of outflows underestimates the number of 
young stars and therefore the star formation efficiency. In 
the L 1495 bowl, a high percentage of young stars will be 
weak-line TTSs that are no longer associated with molec- 
ular or even HH flows. Note also that Taurus in general is 
considered to be a less active region of star formation, being 
associated with less luminous young stars, narrower molec- 
ular line widths and lower kinetic temperatures than e.g. 
Perseus or Orion (Jijina, Myers & Adams 1999). 



4.4 Outflows, cloud structure and the large-scale 
B-field in L 1495 

Taurus has been the subject of many studies of the large- 
scale cloud stru cture and its r elatio nship to the surrounding 
magnetic field; iHeverl l| 19881 ) and iGoodman et~ai1 (|l992l ) 
measure the polarisation of background optical and infrared 
starlight via selective dust absorption, observations which 
illustrate th e field orientati o n in t he ou ter, low column den- 
sity re gions. iTroland et al.l fl996) and ICrutcheT fc Trolandl 
(2000) present complementary Zeeman effect observations 
in 18 cm OH emission. Although these probe higher column 
densities, they only yield information on the line-of-sight 
field component, and are of lower (arcminute) spatial 
resolution. To better trace the B-field that pervades the 
dense cores and high-density gas that envelopes these cores, 
one requires high spatial resolution observations of polarised 
(sub)millimetre continu um emission from magne t ically- 
aligned grains fe.g. [Matthews. Wilson fc Fiegel l200ll; 



Matthews. Fiege fc Moriartv-Schievenl 2002T Crutcher et al.l 
2004! ) . 



From studies of other low- and intermediate-mass 
star forming regions we know that the field orienta- 
tion, projected onto the plane of the sky, often snakes 
through filaments and cores, chang i ng dire c tion on tenths 
of p a rsec scales (e.g. Davis et al. _ 2000bl; Momose et al.l 



120011 ; iMatthews fc Wilson! 120021 ; iHoude et al.1 12002T ). S 

veys of cores and clumps in star forming clouds also re- 
veal a general lack of order in the orientations of the 
long axes of oblate starless and p rotostellar cores (e.g. 
lHatchell et al.ll2005l ; lKirk et al.|[2006l ); similarly, outflow sur- 
veys show that HH and H2 line-emission jets are usually ran- 
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domly orientated ( 


Stanke, McCaushrean & Zinneckerll2002l; 


IWalawender et al.l 


20051; Davis et al.l 20091). Together, these 



data suggest that the magnetic field in the high-density re- 
gions is poorly coupled to the neutrals, allowing ambipolar 
diffusion to build up mass on star forming cores, or that the 
magnetic energy is insufficient to overcome the kinetic en- 
ergy associated with the dense clumps and cores. Either way, 
the orientation of outflows and their associated protostellar 
cores do not appear to be strongly linked to the large-scale 
B-field. 

But could L 1495 be an exception to the rule? The 
south-east ridge in L 1495 is notable for being orien- 
tated roughly perpendicular to the surrounding B-fiel d 
l|Goodman et al.lll992i ; iHeilesI l2000l : iGoldsmith et"afll200Sh . 
Goodman et al. give a mean polarisation position angle 
along the ridge of 27° from optica l data, and 31° from in - 
frared polarisation measurements. IGoldsmith et al. (2008) 
over-plot polarisation vectors (compiled by Hcilcs 2000) onto 
their 13 CO integrated intensity map of Taurus and find that 
the field is not only perpendicular to the L 1495 ridge, but 
is aligned with "striations" in the surrounding low column 
density medium. This suggests that - in these low density 
regions - the B-field is well coupled to the molecular gas 
through ion-neutral collisions (UV penetration maintaining 
a degree of ionisation in the gas). Collapse along these field 
lines may then produce the chain of high density condensa- 
tions that forms the south-east ridge in L 1495. Although the 
field within the dense ridge awaits sensitive submm contin- 
uum polarisation measurements (these are planned as part 
of the Gould Belt survey), we point out here that many of 
the outflows found along the south-east ridge are orientated 
perpendicular to the ridge, and parallel with the field direc- 
tion in the surrounding low column density gas. This obser- 
vation is in general disagreement with the results of Greaves, 
Holland & Ward- Thompson (1997), who in a study of five 
protostars with outflows find that if the flow is in the plane of 
the sky the field tends to be orientated perpendicular to the 
flow axis. However, Greaves et al. used submm continuum 
polarisation measurements and therefore probed the higher 
column density regions close to each outflow source. It will 
be interesting to see whether the field direction, measured 
with the polarimetry facility currently being developed for 
SCUBA-2, changes close to the young outflow sources in 
L 1495. 



5 CONCLUSIONS 

HARP observations in CO 3-2 emission are shown to be ideal 
for tracing outflow activity in nearby star forming regions. 
The observations discussed here reveal as many as 16 molec- 
ular outflows in L 1495; most are associated with HH objects 
and/or molecular hydrogen line-emission features (MHOs). 
Candidate outflow driving sources (protostars or TTSs) are 
identified for eight CO flows, although only four of these 
flow progenitors appear to be associated with H 13 CO + cores. 
Even so, we note that the CO outflow-driving sources have 
redder near-IR colours than their HH jet-driving counter- 
parts. We also find a possible correlation between outflow 
mass and the associated core density and size, the more mas- 
sive flows being associated with the denser, more compact 
cores. 



OMK02 find that 13 of the 22 H 13 CO+ cores in L 1495 
have a column-density mass that exceeds the virial mass. 
Of these 13 cores, five seem to be associated with protostars 
(and four of these drive outflows). We estimate that the ratio 
of prestellar to protostellar cores in L 1495 is approximately 
in the range 1.3-2.3; the ratio of starless to prestellar cores 
is estimated to be ~1. 

Overall, we find that the bowl is more evolved than the 
south-east ridge in L 1495; in the bowl there is a paucity 
of molecular flows though a larger fraction of TTSs. The 
fact that the ridge is spatially more compact (long and nar- 
row) in comparison to the bowl supports this interpretation. 
The star formation efficiency in L 1495 (particularly the 
south-east ridge), estimated from the ambient cloud den- 
sity, a canonical value for the cloud-to-core mass ratio, and 
the observed number of outflows, is consistent with other 
low-mass star forming regions. 

In comparison to Orion and Perseus there is a modest 
number of outflows and protostars in L 1495. However, the 
region, especially the south-east ridge, is relatively simple 
(in comparison to, say, NGC 1333), with little chance for 
source confusion and outflows overlapping other outflows or 
unrelated molecular cores. The statistics, though modest, 
are therefore likely to be more robust than in other regions. 
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